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We show that atomic interference in the reflection from two suitably polarized evanescent waves is 
sensitive to retardation effects in the atom-surface interaction for speciflc experimental parameters. 
We study the limit of short and long atomic de Broglie wavelength. The former case is analyzed 
• in the semiclassical approximation (Landau-Zener model). The latter represents a quantum regime 

I and is analyzed by solving numerically the associated coupled Schrodinger equations. We consider 

^SJ ■ a specific experimental scheme and show the results for rubidium (short wavelength) and the much 

^ ' lighter meta-stable helium atom (long wavelength). The merits of each case are then discussed. 

I. INTRODUCTION 

The interaction between a ground-state atom and a dielectric or conducting wall has been investigated theoretically 
g ' ( i, i, §, i, i, i, i, §, §, 0) and experimentally ( 0, 0, (l|). Theoretical studies have been 
performed on different levels, from a simple model of a dipole-dipole interaction of the atom and its mirror image, to 
the full QED relativistic quantum treatment. Interesting in particular are the long-range Casimir interactions ||^ that 



^ ' were recently observed in cavity QED experiments , . When the atom-wall distance z is not small compared 
, ^ ] to the wavelength of the dominant atomic transitions, the law associated with the instantaneous electrostatic 
^ interaction is no longer valid. The full QED treatment leads to the famous long distance z"^ law. 
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Recent experimental developments enable precise manipulation of cold atoms by lasers, see e.g. Ref. |15[|. Small and 
well defined velocities of the atoms can be achieved using advanced cooling and launching techniques and a detuned 
laser field can be used to create controlled and adjustable potentials for the atoms. Under these conditions, atoms 
can be used to explore the atom-surface potential, for example using evanescent-wave mirrors. Classical reflection 
from such an atomic mirror was used to measure the van der Waals force between a dielectric surface and an atom 
in its ground state . This experiment though, could not fully discriminate between the electrostatic and the QED 
expressions. Segev et al. jl^ considered a similar experiment in the quantum regime (atoms undergoing above- 
I barrier, classically forbidden reflection). Unlike classical reflection, which can only be used to identify thresholds and 
T-H . to measure the height of the potential barrier, quantum probabilities are determined by the complete potential curve, 
and are sensitive to the short- and long-range behavior of the potential. It was found that, for velocities of the order 
, of the recoil speed, the quantum reflection probabilities are indeed sensitive to the long-range (Casimir) interaction. 
■ In this work we study how the form of the atom-surface interaction can be observed using atomic interference of 
the type reported in We consider atoms with multiple ground state sub-levels, which feel different potentials 

. ^ ■ in the evanescent radiation fleld. These potentials can be probed by using stimulated Raman transitions within the 
^ ' evanescent wave. These transitions exhibit interference effects whose phases depend on the atomic trajectories and on 
the entire potential, as in quantum reflection. An important aspect of the effects we discuss here is that they occur 
for higher incident velocities than those considered in and may therefore be experimentally easily accessible. 
Furthermore, even for small velocities (near the recoil velocity) the effects are more dramatic than those shown by 
. ^ ' quantum reflection. 

k> " As discussed in Refs. and jl^, the experiment is analogous to an atom interferometer whose size is of the 
5_j ' order of a fraction of the evanescent wave decay length, We will focus on three experimental cases: in the first, 
' the atomic de Broglie wavelength AdB ^ i^^^ ■ This leads to interference fringes analogous to those in any other 
interferometer. It corresponds to the experiment described in Ref. JTzt in which rubidium (Rb) atoms are dropped 
from a height of about 2 cm. Second, we will repeat the situation of the first case, but with Rb replaced by meta-stable 
helium (He*). In this case we have AdB/27r ^ a few which allows for only a few interference fringes. This allows 
us to illustrate how the signal changes when only the atomic species is changed. Finally, we will examine the case 
where AdB/27r ~ taking He* atoms at the recoil velocity as our example. No fringes are present, but a strong 
dependence on the nature of the atom-wall potential is demonstrated. To simplify the discussion, we work with a 
two-level model, between the initial atomic state and the adjacent one in the ground state manifold, thus neglecting 
the populations of the other ground state sub-levels as well as the excited levels. For the atom-wall interaction we 
will use the published values for each atom. 
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In the short AdB limit, the motion of the atom can be treated semi-classically. We thus calculate the transition 
probability between the two atomic levels according to the Landau-Zener model for adiabatic transitions. In the other 
cases, the study of the atomic motion requires a fully quantum mechanical treatment. We study these cases by solving 
the associated coupled Schrodinger equations numerically. 

In Sec.|| we describe the physical system and the theory for atomic interference. In Sec. Ill the results are presented. 
Discussion and conclusions are given in Sec.LV. 



II. MODEL FOR INTERFERENCE OF MULTIPLE GROUND STATE SUB-LEVELS 

We will focus on an experimental setup along the lines of Fig. |l|. Here, a strong laser beam with frequency lj 
and transverse magnetic (TM) polarization (i.e. magnetic field perpendicular to the plane of incidence) creates an 
evanescent wave, which is nearly a~ circularly polarized with respect to the y-axis. A second weak laser beam with 
frequency lo — Alu and transverse electric (TE) polarization (i.e. electric field perpendicular to the plane of incidence) 
creates another evanescent wave with tt polarization along the y-axis. 

Atoms normally incident on the evanescent wave move in an effective optical potential Might (which in general 
depends on the internal state of the atom) and an attractive atom-wall interaction Vwaii- The total potential is given 
by 

U{z) = VnM^) + V^.ii{z). (1) 



A. Optical potential 

Let us consider an atom in the evanescent field and assume for simplicity that the cr~ polarized wave is strong 
compared to the tt wave. The (t~ component lifts the magnetic (Zeeman) degeneracy of the atomic ground state 
(the quantization axis is in the y direction), so that each magnetic sub- level feels a different optical potential. To 
first order, the tt polarized wave produces a coupling between the atomic sub-levels via a Raman transition: for 
example, starting from the sub-level m^, the atom absorbs a a~ polarized photon and emits a stimulated photon with 
TT polarization. The atom thus ends up in the — 1 sub-state with its total energy increased by hAcu. This transition 
is resonant when hAto is equal to the energy difference between the magnetic sub-levels. If the two evanescent waves 
are counter-propagating, in the reference frame moving with the optical grating, the situation corresponds to grazing 
incidence diffraction [|l7|. For a review of the theoretical understanding of atomic diffraction and interference from 
evanescent waves, see |19[|. 

In the limit of low saturation and a detuning d large compared to both the frequency difference Alu and the 
natural linewidth, the excited state manifold may be eliminated adiabatically, and, for a ground state of total angular 
momentum (including nuclear spin) Jg, the atomic wave function is described by the 2Jg + 1 Zeeman components 
m >,m = — Jg, ...,+Jg. An atom at a distance z from the surface of the mirror is subject to an optical potential 
Might (-z) whose matrix elements are of the form [we suppose that the frequency difference Aw and Zeeman shift are 
negligible compared to the detuning S and the hyperfine structure of the excited level] 

< m\Viif.Y,t{z)\m' >= — ^ E*{z)Eq,{z){Jgm;lq\Jeme){Jeme\Jgm';lq'), (2) 

where a product of Clebsch-Gordon-coefficients appears on the rhs, the electric field polarization is expanded in the 
usual spherical basis with coefficients Eg, [q — —1,0,-1-1) and d is the reduced dipole moment. The optical potential 
couples different Zeeman sub-levels if the field is not in a pure polarization state with respect to this basis, as in the 
setup considered in this work. The optical potential due to the two evanescent waves can then be written as 

(f 

< m\Viightiz)\m' >= —Cjnm' exp(-2Kz), (3) 

where the Cmm' coefficients are given by the Clebsch-Gordon coefficients times the field amplitudes at the surface 
2; = (see Eq. The inverse decay length n \s n = fc[(nsin0)^ — 1)]^^^, where k = uj/c is the free field wave vector, 
assumed to be the same for each laser, n is the refraction index of the prism and 9 is the angle of incidence of the 
lasers with the surface. In our calculations we will always use n — 1.87 and 6 = 53°. The wave- vector k is different 
for each atom. 
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B. Atom— wall potential 



The simplest model for the interaction of a ground state atom and a wall of dielectric constant e considers the 
interaction between a dipole d and its mirror image and yields the Lennard- Jones potential, 

where < (i| > and < > are the expectation values of the squared dipole parallel and perpendicular to the surface 
101 , | |20[ | . This expression for the potential is approximately valid for e independent of frequency and kz much smaller 
than unity. 

If we take into account retardation effects, the Casimir-Polder potential is obtained, where the finite propagation 
time between the dipole and its image results in a different power-law behavior for large z j2| , 

hm VZiliz) cx z-\ (5) 

z — too 

In the complete QED theory the interaction potential between an atom of polarizability a{uj) at a distance z from 
a dielectric wall can be written as H]: 



^QED. N _ 



wall V'J g^2c3 

with 



CarW = -^Z23y^ duJiu^a{iu)x dp + dpjH{p,e)exp{-2tpioz/c) (6) 



H{p) ^ ^ ;-'+P' -P + (1 - 2p^)^^±£^. (7) 
y/e-l+p'^+p y/e - 1 + p'^ + ep 

The numerical values of the constant C^^^ for Rb and He* atoms used in this paper are given in Table |, as well 
as the numerical values of some other important parameters. Since the atom is in a ground state with I = 0, the 
value of C^^^ is the same for any magnetic or hyperfine sub-level. The detailed interaction between a ground state Rb 
or meta-stable He atom and a dielectric surface were recently calculated. For the van der Waals potential we have 
used data from for Rb and from for He*. In the former work an interpolation formula for the van der Waals 
potential is given as V^^f^^^z) = V4fi(z)0.987[(l + 1.098z)-i - 0.00493z(l -I- 0.00987^3 _ 0.00064z'*)-i] where z is 
expressed in units of the laser wavelength A/(27r). This formula approximates the numerical calculation with a 0.6% 
accuracy between and 10A/(27r). 



C. Population transfer 

In this work we will consider only population transfer from the incident sub- level i to the final one / (m/ = — 1), 
that is only two levels. This is a good approximation if the coupling is weak enough for the population of the other 
levels to be negligible. An example of the total interaction potential from Eq. |^ is shown in Fig.y. The Uff potential 
curve has been shifted vertically by —HAuj, corresponding to the kinetic energy change. Then the two (adiabatic) 
potential curves cross at the point of resonance. The coupling turns the exact crossing into an avoided crossing. An 
incoming wave function in the rrii channel is split in two parts that are subsequently reflected from their respective 
repulsive potentials and recombined after the second passage at the crossing. Thus, the evanescent wave realizes a 
"Michelson interferometer" with a single beam splitter and two mirrors. 

When the atomic AdB is short, the avoided crossing can be treated by means of the semi-classical Landau-Zener 
model for non-adiabatic transitions [ p3| . Assuming that the atom moves through the crossing with a constant velocity 
Vc (fixed by energy conservation), the Landau-Zener formula allows one to compute the probability amplitude for the 
two atomic levels after the crossing. The transition probability Wif from the initial sub-level i to the final sub-level / 
is given by 

w^j = ATLzRLzcos'{S(f>), (8) 

where R^z = 1 — T^z and the transmission coefficient is T^z = exp(— ttA), with A = 2| < i\U\f > l"^ / (h"^ kAujVc) . 

The phase difference S(f> is given by the difference in the phase shifts between the crossing and the turning points 
in the semi-classical approximation, plus a correction term [ pif : 
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dzpf{z) - / dzpi{z) 



where Zc is the position of the crossing point, z„ the classical return point for an atom in the n-th level, Pn{z) its 
momentum and T the Gamma function. 

Changing the frequency difference Aa;/(27r) causes a change in the length of one of the interferometer arms, thus 
a change in the phase difference 54> between the two paths. As a consequence we expect the transition probability to 
show oscillations in Aw (Stiickelberg oscillations). We will see in the next section that dtf) is very sensitive to the exact 
shape of the potential. The amplitude of the oscillations also depends on Aw both explicitly and implicitly through 
Vc (the crossing point moves with changing Aw). 

The Landau-Zener model is a good approximation only when the atom speed is approximately constant during the 
crossing. In particular it is not valid when the classical return point and the crossing are close to each other or when 
the de Broglie wavelength of the atom is of the order of the width of the interaction region. In order to explore this 
long-wavelength regime, we have to forgo the semi-classical Landau-Zener model and solve numerically the coupled 
Schrodinger equations for the system. Since atoms that cross the potential barrier stick to the dielectric surface, the 
appropriate boundary conditions at the surface are those for a running wave propagating downward {z — > — oo), while 
for z — > oo the solution is a superposition of downward (incident) and upward (reflected) waves. We have integrated 
the system of Schrodinger equations using the renormalized Numerov method p^ ]. To avoid the singularity at the 
surface we have modified the potential to be a large negative constant near the interface and verified that the transition 
probability does not change by varying the value of the constant. 



III. RESULTS 



We will present our calculations in two parts. First we discuss the short wavelength regime, and point out various 
experimental strategies to observe retardation effects. Then we discuss two cases in which the de Broglie wavelength 
is not small compared to the evanescent wave decay length (the "long wavelength regime"), in which a numerical 
integration of the Schrodinger equation is necessary. 



A. Short wavelength regime 

Fig. H shows a calculation of the population transfer Wif as a function of the frequency difference Aw for the 
Lennard-Jones (LJ) and QED model of the van der Waals interaction. The value of the light potential is the same 
for the two curves (i.e. the laser intensity is the same), and the incident momentum is 115?ifc, which corresponds to 
Rb atoms dropped from a height of 2.3 cm. We use the Landau-Zener approach to calculate the transfer probability 
(see Eq. ||). As in and we observe Stiickelberg oscillations in the transfer probability. These oscillations 
can be understood as the variations in the accumulated phase difference between the two different paths taken by the 
atoms after the level crossing shown in Fig. 2. The de Broglie wavelength AdB is such that several fringes appear as 
the position of the level crossing is moved through its possible range. The last oscillation at the higher frequencies 
is where the crossing point and the classical return point are very close to each other and the Landau-Zener model 
breaks down. We have set the probability to zero beyond this limit. In reality, the transition probability falls roughly 
exponentially with frequency, as one finds solving the Schrodinger equations numerically (see next subsection). 

In general, we find that the dephasing between the two curves (with and without retardation effects) is greatest 
when the atoms are incident at an energy close to the top of the potential barrier. Note that the height of the barrier 
is greater when retardation effects are included, since these reduce the strength of the atom-surface interaction. The 
effect of retardation is roughly to shift these fringes by half a fringe to the left. The major cause of this shift is 
the increase in the height of the total potential which is greater for the i level than the / level. A 10% increase in 
the value of the light-shift potential would exhibit nearly the same shift. Therefore, since it is not possible to turn 
retardation on and off, it would be necessary to measure absolutely the light-shift potential to better than 10% in 
order to distinguish a retardation effect. Experimentally this is rather difficult. Instead of attempting to measure 
the absolute light shift, one could rather measure the absolute height of the potential by observing the threshold of 
reflection as in ||ll[, and using the known kinetic energy of the atoms to get an absolute calibration of the height. Let 
us assume then that the barrier height, instead of the light intensity, is known. In Fig.^a we show the result with the 
same parameters as in Fig.^ except for the light intensity in the LJ model, which has been changed so as to have the 
same barrier height as the QED model. In this case the shift is much smaller, about 1/5 of a fringe. We have verified 
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that even taking into account an experimental uncertainty of a few % in the height of the corresponding potentials, 
the two models are still clearly resolved. 

This approach seems feasible, but a third method of observing the effects of retardation is possible if one uses more 
of the information available in the oscillation curve. Fig. 4b shows the same curves as Fig. 4a but with the QED 
curve numerically shifted so as to coincide with the LJ curve at its maximum. One sees that the period of these 
oscillations is not the same. It decreases with decreasing detuning, faster for the full QED potential, so that there is 
a difference in the spacing of the minima in the population transfer. Thus with fringes with sufficiently high signal 
to noise, one can distinguish retardation while leaving the absolute barrier height as a free parameter in a fit to the 
data. It seems to us that a viable experimental method is to use a combination of the second two approaches. Careful 
measurements of the barrier height can be used to cross check a fit to the Stiickelberg oscillations with the barrier 
height as a free parameter. 

The incident energy, momentum and barrier height used in Figs. 3 and 4a were arbitrarily chosen to correspond to 
the experiment in p7[ |, but it would be interesting to know how Figs. 3 and 4a would change, especially for different 
incident momenta (de Broghe wavelengths). We studied this question by repeating the calculation of Fig. 4a, for 
different incident momenta, while always keeping the barrier height 10% above the incident kinetic energy. We find 
that, roughly speaking, the number of oscillations increases as the incident momentum. This is because the number 
of fringes in the interferogram increases with decreasing wavelength. The accumulated phase difference between the 
oscillations in the LJ and QED models, over the corresponding frequency range, also increases approximately linearly 
with incident momentum. Thus if we consider the fringe shift divided by the fringe period as a figure of merit, the 
sensitivity of the experiment to retardation effects increases with increasing incident momentum. 

B. Long wavelength regime 

We now examine the large de Broglie wavelength limit, where the semi-classical model breaks down since the atomic 
wavelength is not small compared to the interferometer size. We first consider the case of He* atoms dropped from 
a height of a few cm as for Rb. We show the results in Fig. ^. For an initial distance of 2.3 cm from the mirror, 
the incident momentum of He* is 7 Ahk, much lower than for Rb. Again we will chose the intensities of the strong 
a~ wave so as to have the same barrier height for the two potentials, about 10% above the incident kinetic energy. 
In this regime, to observe interference, the detuning between the evanescent waves must be of only a few MHz, in 
fact, beyond about lOMHz the crossing point is closer to the surface than the return point. One only sees one or 
two Stiickelberg oscillations since the momentum involved is small and the atoms do not accumulate enough phase 
difference to show more oscillations. The two potentials give similar results (Fig.^, the main difference being in the 
shape and height of the big peak. 

We have also looked at even lower incident momenta. Near the recoil limit, one can expect interference only for 
detunings of less than one MHz and the accumulated phase difference is too small in this range to show any Stiickelberg 
oscillations. Both for He* and Rb one obtains a big peak in the transition probability and no Stiickelberg oscillations, 
as expected. The difference between the van der Waals potential and the full QED potential still shows up in the 
different shape and height of the peaks (Fig.^) . Here the width of the curve is delimited by the frequency for which 
the classical return point and the crossing point coincide. 

Finally we note that the qualitative features of the short and long Ads limit are the same for Rb and He*, e.g. He* 
with an incident momentum of about 100?ifc gives the same type of oscillation pattern as Rb. 

IV. CONCLUSIONS 

In summary, we have proposed an experiment to probe van der Waals like surface interactions by exploiting 
interference mechanisms for well defined Zeeman sub-levels of atoms moving in two evanescent waves. Retardation 
can be resolved using atoms incident at speeds which are easily obtained in free fall over a few centimeters. The 
controlling parameter is here the detuning between the two evanescent waves. One then measures the fraction of the 
atoms which have undergone a change of magnetic sub-level as a function of this detuning. 

For the situation in which the atomic de Broglie wavelength is sufficiently small the experiment resembles typical 
interferometry experiments. The theoretical description is semi-classical, employing well defined atomic trajectories, 
while experimentally, one seeks a particular (non sinusoidal) fringe pattern as a signature of retardation. This should 
be possible with an improved version of the experiment of Ref. . 

Another approach is to investigate the interaction of atoms whose de Broglie wavelength is not small compared 
to the length scale of the interferometer. In this regime most of the information is to be found in the shape of the 
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population transfer curve, since there are very few or zero interference fringes. Note though, that we have assumed 
throughout that the incident atoms are mono-energetic. This means that the velocity spread of the incident atoms 
must be small compared to the atomic recoil. Nevertheless this may be worth the effort, because the predicted effect 
of retardation is quite dramatic. 

For a quantitative comparison however, the presence of all the sub-levels (which give rise to the multiple crossing 
of the dressed potentials at the same distance z,, from the surface) and, possibly, losses from spontaneous emission 
have to be taken into account, but the results are not qualitatively different. 
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Parameter Rb He* 



A 


780 X lO^^m 


1083 X lO^^m 


r 


2tv5.9 X 10*^8-1 


27rl.6 X lO'^s-^ 




6.4 X iQ-^nr 


2.6 X io-2;ir 




o.iis^ir/A^ 


0.125?ir/A^ 


state i 


5Sy2,F = 2,mF = 2 


2^Si,J = l,mj = 1 


state f 


551/2, -F = 2, mp = 1 


2^Si,J = l,mj = G 


Cii 


\E-i{0)\H/3 
\E.i{0)\H/2 


|iJ-i(0)pl/6 
\E-m?l/2 


Cff 


Cif 


E^^{OyEo{0)^/2/3 


£_i(0)*i;o(0)l/3 



TABLE I. Values of parameters used in the text for Rb and He*: laser wavelength A, atomic transition width F, recoil 
energy Er, van der Waals coefRcient C''^' (see Eq.^). Definition of the initial and final atomic state (i and /). Optical potential 
coefficients dj (see Eqj^. We assume here that the cr" wave (-B-i) is much stronger than the tt wave {Eq). The dielectric 
constant of the wall is e = 3.49 and n — 1.52fc. 
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Atoms before 
the bounce 



Probe beam 
— > 



Atoms after 
the bounce 




FIG. 1. Typical experimental setup. The atoms are released from a MOT above the prism. Two slightly detuned laser beams 
of polarizations TM (and frequency uj/{2n)) and TE (and frequency {lu — Auj)/(2ti)) form evanescent waves on the surface of 
the prism of polarizations a~ and n respectively (y is the quantization axis). The reflected atoms are detected by a detection 
beam. 



800 




FIG. 2. Typical potential curves experienced by the atoms during reflection. The potential for state / is shifted by the 
Raman energy —hAcu. The atoms approach the potential barrier in state i, pass through the curve crossing twice and can end 
up on either state i or /. For each final situation, two paths are possible (between the crossing and the turning points) and 
can interfere producing fringes as a function of the location of the crossing. 
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5 10 15 20 25 30 

Frequency difference (MHz) 

FIG. 3. Population transfer from state i to state / after reflection (wif from Eq. ^) vs frequency difference Alj/{2tv) for 
rubidium atoms, released 2.3 cm above the mirror, (i.e. with incident momentum llSRfc). The solid line is for the QED model 
and the dashed line for the Lennard-Jones (LJ) model. The coupling coefflcients defined in Eq. |^ are \Cif\ — S59E^^ and 
\C,,\ = 3.8 X W^E^^ 
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FIG. 4. Upper figure: population transfer Wif (see Eq. |8|) vs frequency difference Atj/(27r) for rubidium atoms released 2.3 
cm above the mirror (i.e. with incident momentum 115?ifcRb). The solid line is for the QED model and the dashed line for the 
Leenard-Jones (LJ) model. The height of the potential barrier is 1.48 x 10'' -E^'^, i.e. the intensity of the strong laser beam has 
been adjusted to give the same barrier height for both models. The light-shift coefficients for the QED model are the same as 
for Fig. |, while for the LJ model |Cii| = 4.11 x W^E^^ and \C\f\ = 890E^^. Lower figure: same parameters as in the previous 
figure with the LJ curve artificially shifted 0.79 MHz to the left in order to show the changes in fringe spacing. 
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FIG. 5. Population transfer Wif (see Eq. ^ vs frequency difference Alu/{2iy) for meta-stable fielium atoms released 2.3 cm 
above the mirror (i.e. with incident momentum 7.4hkuc)- The solid line is for the QED model and the dashed line for the 
Lennard-Jones (LJ) model. The height of the potential barrier is Q1E^° , i.e. the intensity of the strong laser beam has been 
adjusted to give the same barrier height for both models ( 10% above the incident kinetic energy). The light -shift coefficients 
for the QED model are \Cii\ = 252£?° \Cif \ = 20.8£?'= while for the LJ model \C^i\ = 303£"° and \Cif \ = 22.7£^'°. 
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FIG. 6. Population transfer Wif vs frequency difference Alj/{2-k) for meta-stable helium atom with incident momentum 
hkiie- The solid line is for the QED model and the dashed line for the Lennard-Jones (LJ) model. The height of the potential 
barrier is 1.28£'^°, i.e. the intensity of the strong laser beam has been adjusted to give the same barrier height for both models. 
The light-shift coefficients for the QED model are \Cii\ = O.SE^" = OAQE^" while for the LJ model \Cu\ = 1.47S^° and 



11 



